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Abstract 

J/ip production cross section considering the contribu- 
tions of both color-singlet and color-octet cc channels is 
calculated. The result is used to study the suppression 
of J/ip in nuclear collisions. With absorption cross sec- 
tions for (cc)8 (Jabs — 11"^^ and (cc)i cr^^^ ~ Om6 the p-A 
and A-B data except for Pb-Pb can be explained. Possi- 
ble explanations of additional suppression in Pb-Pb are 
discussed. 
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I. INTRODUCTION 



It was expected that a suppression of J/ip production in relativistic heavy ion 
colhsions can serve as a clear signature for the formation of a new matter phase 
Quark Gluon Plasma (QGP) [jl|. This suppression effect was observed by NA38 col- 
laboration later 0]. However, successive research pointed out that such suppression 
could also exist in hadronic matter (HM), even though by a completely different 
mechanism 0]. The anomalous J /ip suppression recently reported by the NA50 col- 
laboration and there have been a number of attempts to explain it, such as the 
onset of deconfinement, hadronic co- mover absorption and the energy loss model 
[|^. To understand the expermental data clearly, the formation and absorption 
mechanism of J/ip must be studied carefully. 

Quarkonium production has traditionally been calculated in the color singlet 
model. However, it has become clear now that the color singlet model fails to provide 
a theoretically and phenomenologically explanations of all production processes and 
the inclusion of color octet production channels removes large discrepancies between 
experiment and the predictions of the color singlet model for the total production 
cross section. In principle, the J/ ip state is described in a Fock state decomposition 

|J/^ > = 0(l)|cc(=^5f^) > +0{v)\cc{^Pf)g > +0{v^)\ccCsll'^)g > + 

0{v')\ccCs[''%g > +0{v')\ccCD^j'%g > +..., (1) 

where ^s+i^O^.s) characterizes the quantum state of the cc with color-singlet or octet 
respectively. This expression is valid for the non-relativistic QCD (NRQCD) frame- 
work and the coefficients of each component depend on the relative three-velocity 
1^1 of the heavy quark, and under the limit of \v\ 0, i.e. both c and c remain 
relative rest, eq.(|I|) recovers the expression for color-singlet picture of J/ip where 
0(1)^1. 

The color-octet component {cc)s seems to play an important role in interpreting 
the Collider Detector at Fermilab (CDF ) experimental data [0], and further studies 
show that it may also have great influence on the quarkonium production at other 
collider facilities P|-pi^. The investigations [|T^ on quarkonium hadrproduction at 
fixed target energies find that color-octet contribution to the cross section are very 
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important, and the inclusion of color-octet production channels removes large dis- 
crepancies between experiment and the predictions of the color-singlet model for the 
total production cross section. Therefore one can expect that the color-octet {cc)s 



would also manifest itself in heavy ion collisions [|T5|. From the above discussion, 
one can know that the charmonium production can be devided into two steps. The 
first step is the production of a cc pair. The cc pairs can be either (cc)i or {cc)s. The 
cc pairs are produced perturbatively and almost instantaneously, with a formation 
time Tf ^ (2mc)^^ ~ 0.07 fm in the cc rest frame. The second step is the formation 
of a physical states of J/ip, that need much longer time. People believed now that 
J/ ip suppression in hadron matter can be consider as pre-resonance absorption. Satz 
first proposed to use the pre-resonance absorption model to explain nuclear colli- 



sions data |15|. However, in their work the pre-resonance state of charmonia is only 
in color octet. In this paper we give a complete leading order calculation of charmo- 
nium production cross section for (cc)i and (cc)8 in corresponding energies. Then 
the results is used to study charmonium suppression in heavy ion collisions. As it is 
known in the hadron-hadron collisions, the (cc)i produced is almost point like, and 
the color-octet can interact with nucleons much more strongly than the color-singlet 
(cc)i, so the color-octet would dissolve much faster into D and D than (cc)i, i.e. 
^abs ~ << (7^5^. Thus the absorption of (cc)i and (cc)8 should be considered differ- 
ently. We find that in our model with the absorption cross section ali^^ ~ Omb and 
a^fjg ~ limb the p-A and A-B data except for Pb-Pb can be explained. To explain 
Pb-Pb data, additional sources causing nonlinear suppression in nuclear colhsions 
must be onset. 

Our paper is devided into four parts. In section II, we describe the color-octet 
scenario borrowed from pp collisions and calculate the charmonium production cross 
section for (cc)i and (cc)8. In section III, the pre-resonance nuclear absorption model 
considering both (cc)i and (cc)8 are described. In the last section the result of our 
model is compared to the expermental data. Some possible models for explaining 
the anomalous J/ip suppression in Pb-Pb collisions of NA50 experiment, such as the 
formation of QGP are discussed. 



3 



II. FORMULATION 



According to the NRQCD factorization formalism, the inclusive production rate 
of heavy quarkonium H in parton level can be factorized as 

a{ij ^H) = Y. a{ij ^ QQ[n]) < O^H > . (2) 

n 

Here, a{QQ[n] + X) describes the short distance production of a QQ pair in the 
color, spin and angular momentum state n, which can be calculated perturbatively 
using Feynman diagram methods. < >, the vaccum expectation value of a 

four fermion operator in NRQCD ||^, describes the nonperturbative transition of the 
QQ pair hadronizating into the quarkonium state H. The relative importance of 
the various contributions of n in eq.(0) can be estimated by using NRQCD velocity 
scaling rules. An important feature of this equation is that QQ pairs in a color-octet 
state are allowed to contribute to the production of a color singlet quarkonium state 
H via nonperturbative emission of soft gluons. Accordingly, the production cross 
section for a quarkonium state H in the hadron process 

A + B H + X (3) 

can be written as 

(Th = Y1 j dXidX2fi/A{Xi)fj/B{x2)(T{ij H), (4) 

where the parton scattering cross section is convoluted with parton distribution 
functions fi/A and fj/B, and the sum runs over all partons in the colliding hadrons. 

At leading twist and at leading order in a^, the color-singlet QQ production 
subproceses for ^'^+^Lj state are 

gg ^ 'So, 'Po,2, (5) 
gg^ 'Si+g, 'Pi + g, (6) 
gq^ ^P, + q, (7) 
qq^ ''Pi+g. (8) 

The corresponding formulae of the above precesses can be found in refs. [0 and 
Tsj . Because the radiative decays Xi,2 J/ip + •j are known to have a large 



branching ratios to J/%1) and the feeddown of the ■0' to J/%1) is also important, their 
contributions should be included in the calculation in reproducing the fixed target 
expermental data of prompt J/%1) production. 

Under the NRQCD factorization scheme, to calculate the quarkonium production 
cross section, one use a double expansions: the perturbative expansion of the short 
distance production amplitude in strong coupling constant and the expansion 
of the nonperturbative long distance hadronization amplitude in typical velocity 
of heavy quark inside the heavy quarkonium. At leading order in perturbative 
theory and up to next-to-leading order in the velocity expansion, the subprocesses 
for leading-twist J/%1) production through color-octet intermediate states are 



qq ^ cc[8,3 S^] ^ J/^ + X, (9) 

gg^cc%^S^]^J/ilj + X, (10) 

gg^cc%^Pj\^Jl'^ + X, (11) 

qq ^ cc[8,=^ Pj] ^ Xj + ^ ^ ^/^ + 7 + X. (12) 

Their cross sections are proportional to the NRQCD matrix elements 

<Q\Oi'^{^Si)\Q>^my, (13) 

<0|C»^/^%)|0>~my, (14) 

<Q\Oi''^{^Pj)\^>^mlv\ (15) 

< 0|Of (2^1)10 >~ (16) 



It is obvious that the above matrix elements are higher order in compared to the 
leading color-singlet ones, but their corresponding short-distance processes are lower 
order in than that in color-singlet processes, this makes the color-octet processes 
non-negligable in reproducing the fixed target experiment data. 

For ip' production the cross section does not receive contributions from radiative 
decays of higher charmonium states, the a^i differs from the direct J/ ip production 
cross section (T{J/il^)dir only in the replacement of V^' matrix elements by J/'0 matrix 
elements. 

With the understanding of quarkoium production, in the calculation we can 
assume that the color-singlet and octet cc are produced at early stage of heavy ion 
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collisions and later evolve according to the environment. Here we assume that the 
cc are produced at the early stage of the ion collisions through hard N-N collisions 
when the phase transition has not yet occurred, and then the system would either 
remains in HM or a new phase QGP is formed at later time. Thus we only need to 
deal with the evolution of (cc), while the production of cc by the deconfined quarks 
or gluons in QGP is ignored. 



III. ABSORPTION IN HADRONIC MATTER 

The conventional probability for & J/ip produced in a p-A collision is given by: 

1 CFpA 



A app 



(fbdzpAib, z)expl^-{A - 1) ^ dz pA{h, z')aabs^ 
= exp{-LApocrabs) (17) 

where app and apA are the J/ip production cross section in proton-proton collisions 
and proton-nucleus collisions, respectively, pA is the nuclear density distribution. 
aabs is an absorption cross section. L is the effective length of the J/ip trajectory. 
It can be derived as 

L = - — — — tqA^^^, for heavy nucleus 

1 A — 1 

= — tqA^^^-^, for light nucleus, (18) 

2 A Tq 

where po = 0.14/m~^ and tq = 1.2/m, Tq = 1.05/m 

As cc pairs are produced almost instantaneously and the formation of a physical 
states J/ip need a much longer time, people now believe that J/ip suppression in 
hadronic matter can be considered as an absorption of pre-resonance cc pairs. As 
discussed in former section, there are both (cc)i and (cc)8 pairs. The (cc)i produced 
is almost point like. The color-octet can interact with gluons much more strongly 
than the color-singlet (cc)i, and therefore would dissolve much faster into D and 
D than (cc)i. Thus their absorption cross section are different. Considering these 
facts, we extend the Eq.([T7|) to 
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p 1 ^PA 

Ka — 



A app 



= /i y (i^bdzpAib, z)exp^ — {A — 1) J dz'f 
+ fsj d^dzpAib, z)expl^-{A - 1) y" dz'pAib, z')al,^'j (19) 

Where fi, fs are relative fractions of (cc)i and {cc)s. o'lbsj'^abs ^^e absorption 
cross sections for (cc) i-nucleon and (cc)8-nucleon correspondingly. 



IV. RESULTS AND DISCUSSIONS 

Before embarking on the computation of cross sections the parameters used in the 
computation should be fixed up. The value of < OlCg''''' (^5*1)10 > can be obtained by 
fitting the theoretical predictions to the CDF Collaboration data at large pt- The 
number of independent matrix elements can be reduced by using the spin symmetry 
relations up to corrections of order v'^ 

< 0|C>f (3Pj)|0 >= (2J+ 1) < 0|C»f (^Po)|0 >, (20) 
< 0\Oi^^CPj)\0 >= (2J+ 1) < 0|08^/'^(=^Po)|0 >, (21) 

< 0|Cf (=^5i)|0 >= (2J+1) < 0\Of'(^Si)\0 > . (22) 

Therefore, the matrix elements < 0|(9|^(^S'o)|0 > and < 0\Og{^Po)\0 > enter fixed 
target production of J/ and ip' in the combination 

As{H) =< 0\OiCSo)\0 > +\ < 0|08^(=^Po)|0 > • (23) 

Up to corrections in f ^, the relevant color-singlet production matrix elements are 
related to radial wave functions at the origin or their derivatives, 

< 0|Of (3^i)|0 >= ^\Rm\ < 0|Of ('Po)|0 >= 7^|P'(0)|^ (24) 

which can be determined from potential model or from quarkonium leptonic decays. 
The values of these parameters, which we use, are ||12| 



7 



< 010^^^(351)10 >= 1.16 G'eV^^ <0\Os^^CSi)\0>=6.6xlO-^ GeV^ (25) 

<0\Ot'CSi)\0>=0.76 GeV^ < 0\Ot' CSi)\0 >= A.6 x 10''^ GeV^, (26) 

< 0|Of (3Po)|0 > /m^ = 4.4 X 10"^ Gey^ < 0|Cf (^^OlO >= 3.2 x 10"=^ GeV\ (27) 

A8(J/^/') = 3.0 X 10"^ GeV^, As{ip') = 5.2 x 10"^ GeV\ (28) 

We use the Gliick-Reya- Vogt (GRV) leading order (LO ) |TB| parameterization for 
the parton distributions of the protons. The c quark mass is fixed to be = 1.5 
GeV and the strong couphng is evaluated at the scale /i = 2mc, that is as ~ 0.26. 
The results of the intergrated cross sections for color singlet and color octet channels 
at several different energies are listed in Table I. 

In the equation (|^) there are two parameters cr^j,^ and af^^, which is different 
from Satz's model |jl5|. As cci produced is almost point like, cr^^^ is very small. 
Thus we take a^^^ = and the value of af^^ is considered as an open parameter and 
determined such as to get the best agreement with the data. In Fig.l we see that 
with (T^jg = llm6 we get quite good agreement with the data except for Pb-Pb data. 
So our pre-resonance absorption model of color singlet and color octet can explains 
both the proton-nucleus and nucleus-nucleus data up to S-U system. To explain the 
non-linear suppression of J/ip additional source of charmonium in nuclear collisions 
must be included. 

The anomalous suppression has been interpreted as a hint of QGP |^. If the 
QGP is formed, the interaction region consists of a hot part where a J/tp can be 
dissociated. Ref |^ modelled the effect of quark-gluon plasma formation by assuming 
that the J ftp produced in those region is completely destroyed whereever the density 
of the participants exceeds a critical value. The survival probability then becomes 

Rab 

= fi J d%(fbAdzAd%BdzBPAibA, ZA)t(h - bA - bB)pB(&B, ZB)9{nc - Upih, bA, be)) 

exp|-(A - 1) J^^ dz'^pAibA, z'^)(jlf,^Y^p[-(B - 1) dz'j^PB{bB, 4)^^afes} 

+/8 J d%d%AdzAt{h - bA - h^)pA{bA, ZA)d%BdzBPB{bB, ZB)9{nc - np(b, bA, be)) 
exp|-(A - 1) dz'APA{bA, 2;^)cT^b,|exp|-(5 - 1) j^^ dz'^PBibB, 4)f^L}> (29) 
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where t(b — ba — be) is called the thickness function, is the density of participants 
per unit transverse area, ric is the critical value. If rtc is chosen to be the highest 
value attained in S-U collisions, i.e., 3.3/m~^ as suggested in Ref P, only in 66% 
part of the interaction region the J /if) can survive [^]. Using this result, in Fig. 2 we 
give a rough analysis of our model with QGP and it is compared with the result of 
our nuclear absorption model. From Fig.2 we can see that using the mechanism of 
QGP, the anomalous suppression in Pb-Pb can be explained. 

However there have been other alternative explanations of the Pb-Pb data by 
NA50, such as collisions with secondary hadrons (co-mover) and the energy loss 
model [0]. To distinguish these mechanisms much more work should be done in 
this field. 
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Figure Captions 

Fig.l, The J/ip survival probability obtained using Eq.(|19]) with crfj,^ = limb 
and a\^g = Omb compared to the expermental data of p-A and A-A collisions at 
different energies. 

Fig. 2, The J/ip survival probability obtained using our nuclear absorption model 
for p-A and A-A collisions at E = 200GeV and the result of Pb-Pb E = 158GeV 
without and with QGP formation in 34% part of the interaction region compared 
to the expermental data. 

Table Caption 

Table I. The intergrated cross sections for color singlet and color octet processes. 
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Table I 



E=450 GeV E=200 GeV E=158 GeV 
(Ji 49.64 nb 24.56 nb 19.54 nb 
(78 95.66 nb 55.48 nb 45.94 nb 
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